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This article studies Markovian stochastic motion of a particle on a graph with finite number of 
nodes and periodically time-dependent transition rates that satisfy the detailed balance condition 
at any time. We show that under general conditions, the currents in the system on average become 
quantized or fractionally quantized for adiabatic driving at sufficiently low temperature. We develop 
the quantitative theory of this quantization and interpret it in terms of topological invariants. By 
implementing the celebrated Kirchhoff theorem we derive a general and explicit formula for the 
average generated current that plays a role of an efficient tool for treating the current quantization 
effects. 
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I. INTRODUCTION 

This is the first in the series of two articles devoted 
to the properties of classical mesoscopic stochastic sys- 
tems under the influence of adiabatically slow driving. 
At thermodynamic equilibrium, laws of thermodynam- 
ics do not allow such systems to perform a directed mo- 
tion on average and only random thermal fluctuations 
can happen. Periodic changes of parameters, however, 
induce pump and ratchet effects [1] that lead to a di- 
rected motion of a system. It is often the case that the 
average currents in such systems are quantized with the 
quantum numbers being integers or fractional numbers. 
This phenomenon was initially discussed for applications 
to mesoscopic quantum mechanical systems at zero tem- 
perature [2-7]. More recently, quantization of currents 
in purely classical mesoscopic systems has also attracted 
attention. It has been discussed for a number of applica- 
tions, including electronic turnstiles [8, 9], ratchets [10], 
molecular motors [11-16], and heat pumps [17]. 

Recently it has been proposed [15] that this phe- 
nomenon occurs generically in stochastic kinetics, where 
we have argued that Markovian evolution of a particle 
in an arbitrary graph leads either to integer or fractional 
quantization of the average current, generated over long 
time in the limit of adiabatically slow parameter driving, 
followed by the low-temperature limit (the latter limit 
is also equivalent to the condition that the amplitude 



of perturbation has much larger energy than ksT). In 
the first (present) manuscript of the series we consider 
in some detail the quantization of average currents phe- 
nomenon, whereas the second manuscript focuses on the 
full counting statistics of currents in the regime of aver- 
age current quantization, and relate the latter to specific 
topological invariants. 

There are numerous reasons to explore the (fractional) 
quantization of currents in stochastic systems. Models of 
periodically driven stochastic systems have applications 
to chemical kinetics and electronics [1]. They have be- 
came a major research direction in the field of stochas- 
tic thermodynamics, which explores the performance of 
microscopic engines beyond standard macroscopic ther- 
modynamics framework [18] . Quantization of mesoscopic 
systems response to periodic driving provides an oppor- 
tunity for robust control of these systems, e.g. for cre- 
ating a source of well-controlled current [8]. In second 
article of the series, we show that the quantization of 
currents is related to some fundamental topological prop- 
erties of the full counting statistics (FCS) of generated 
currents. FCS is an active field of research in mesoscopic 
physics with applications to single molecule experiments 
and nanoscale electric circuits [3, 19]. FCS measure- 
ments can reveal information that cannot be obtained 
by looking at average characteristics. For example, it 
is known that full counting statistics can be related to 
mutual quantum information [20] and can show unusual 
phase transitions in response to periodic driving of pa- 
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ranictcrs [21-23]. Exact results for currents and current 
fluctuations in systems driven periodically by arbitrarily 
strong and fast driving protocols have been established 
[24-30] . Some of the exotic phenomena, such as charge 
fractionalization, can appear in the counting statistics of 
the most simple stochastic models [11, 15, 31]. We shall 
identify and study a topologically nontrivial structure in 
the full counting statistics. Topological phenomena in 
physics arc a subject of a ongoing interest. Studies of av- 
erage quantized currents can result in the classification 
of such new topological phenomena in statistical physics. 

We will also show that the master operators that de- 
scribe the evolution of currents and populations in our 
models are Hermitian up to a transformation that does 
not influence the topological properties of their eigen- 
states. This provides a bridge between the stochastic 
phenomena appearing in our work and quantum me- 
chanics. For example, it implies that all topologically 
protected numbers that we flnd in stochastic models, in 
principle, can be found in some quantum mechanical sys- 
tems. Classification of quantum mechanical pump effects 
is currently an important research direction [32] and our 
work provides a classification for a subclass of such sys- 
tems. 

The main idea that stands behind the mechanism of 
integer current quantization in the low-tcmpcraturc adia- 
batic limit can be described qualitatively in the following 
way. When the time-dependent rates kia, of a transi- 
tion from node i via link a of our network axe in de- 
tailed balance at any time, they can be parameterized 
[see Eq. (4)] by the node energies Ei and barrier heights 
Wa, associated with the network nodes and links, respec- 
tively. In the case of adiabatic driving at low tempera- 
tures the particle sits at the node i with the lowest value 
of Ei, provided the energy is non-degenerate, and there- 
fore no current is generated while the lowest-energy state 
is non-dcgcncrate. Stated differently the currents can be 
generated only when the lowest energy Ei gets replaced 
by Ej for some other node j, which necessarily involves 
passing through a degeneracy Ei = Ej . Passing through 
such a degeneracy along a driving protocol is typical, 
since the latter is represented by a closed curve in the 
parameter space. The current, responsible for the pop- 
ulation flow from node i to node j, is generated around 
the Ei = Ej degeneracy. However, if the barrier val- 
ues arc non-degenerate during the aforementioned pop- 
ulation redistribution, in the low-temperature limit the 
particle will avoid the higher barriers, so that its motion 
will be restricted to a particular spanning tree of our net- 
work, obtained by withdrawing the links with the highest 
barrier values, as described in subsection VIA. There- 
fore, any path that connects the node i to j is restricted 
to the tree. Any such a path can bounce many times 
through links on the tree but contributions to currents 
from transitions back and forward through the same link 
cancel. Hence the contributions to current, produced by 
paths connecting i to j on the tree, arc the same. We 
will say that such paths, which are different only by mul- 



tiple transitions through some links in both directions, 
belong to the same topological class. Since simultaneous 
degeneracy of the node energies and barrier heights is 
non- typical (by the dimensional arguments), for a generic 
driving protocol the population redistribution involves 
the paths of the same topological class. Therefore, when 
over a period of the driving protocol the particle returns 
to the starting node, the current over the period is gener- 
ated by closed paths of the same topological class. Since 
the generated current depends only on the topological 
class of the corresponding stochastic path, all relevant 
stochastic paths, i.e., the ones that provide a substan- 
tial contribution to the current can be obtained by one 
representative. Stated differently, the integral current Q 
(i.e. the flux) that is generated during a period of driv- 
ing may be viewed as having been produced by a single 
closed path and is therefore integer valued. 

To handle a more involved and technical case of frac- 
tional current quantization, which corresponds to a sit- 
uation with permanent degeneracy of some energies and 
barriers, we will map the problem of identification of the 
average instant current J (so that the flux Q is given 
by the integral of J{t) over the period of driving) onto 
the problem of finding the currents flowing over the links 
in an electrical network, given the values of the external 
currents associated with the nodes; the described prob- 
lem is known as the Kirchhoff problem [33, 34]. This will 
allow us to apply the Kirchhoff theorem [33, 34] that pro- 
vides an explicit solution of the Kirchhoff problem. We 
will show below that this solution is a powerful and con- 
venient tool in studying the low-temperature adiabatic 
limit. 

The manuscript is organized as follows. In section II 
we introduce the key objects that play an important role 

in both papers of the series. Section 111 contains a pre- 
liminary discussion of the main results of the present 
manuscript for arbitrary connected finite graph. In sec- 
tion IV, we develop an analogy between the stochastic 
pump effect and currents in an electrical network. Before 
we proceed to the derivation of the main results for an 
arbitrary graph, section V works out in a self-consistent 
way a simple "toy-model" that we use for illustrating 
the basic concepts and phenomena. Sections VI and VII 
provide the explicit derivation of our main results: ex- 
pressions for the integer and the fractionally quantized 
currents, respectively. Section VIII is devoted to topo- 
logical interpretation of the (fractionally) quantized cur- 
rents in terms of the winding index around elements of 
what we call the space of bad parameters. We summarize 
our findings in section IX. 



II. KEY OBJECTS AND TERMINOLOGY 

In this section we introduce the key objects appearing 
in both manuscripts of the series. 
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FIG. 1. A closed graph with 5 nodes and 7 hnks, representing 
a Markov chain with five kia = fce^'^'~'^°'. Nodes are char- 
acterized by energies of weU Links are characterized by sizes 
of barriers Wa, where a = 1, . . . ,7. 

A. Markov-chain model for stochastic pumping 

We start with describing a model of stochastic driv- 
ing, which basically follows the instantaneously detailed- 
balanced Markov-chain model described in [30] with some 
slight differences. In this instance, the temperature will 
be treated explicitly. We will also allow more than one 
link that connects two nodes, as well as links that connect 
a node to itself (self-links). More formally, our network 
is described by a non-oriented graph X = (Xq^Xi) with 
Xq and Xi being the sets of nodes and links, respectively. 
The non-oriented boundary map, d, associates with a link 
a its non-oriented boundary da — {i, j} (or da = {j} in 
the case of a self-link), represented by the set of nodes, 
the link is attached to, as we show in Fig. 2. It is very 
convenient to fix some arbitrary orientation on the graph, 
i.e., put an arrow (direction) on each link. This allows 
oriented boundaries of links to be introduced, which we 
will represent by standard brackets, i.e. da — {doa, c?ia), 
so that if a is not a self-link, the arrow on it goes from 
doa to dia] for a self-link doa — dia. 

A stochastic process is described by a set kia of transi- 
tion rates. The rates can be (and in the case of stochastic 
driving actually are) time-dependent. Stochastic dynam- 
ics can be described by the master equation for distribu- 
tions p = {pj\j e Xq) that has the form 

dtpit) = H{t)p{t), dt = d/dt, (1) 

where the master operator H is given by 

{Hp)i = ^ ^ kjaPj - hpi, (2) 
with the departure rate 

da={i,j} 

ki = ^ ^ ^ ^ kia- (3) 



Here the notation, X^asxi*'"''^ should be understood as 
the sum over all links having the boundary represented 
by nodes {i} and {j}. 

B. Detailed balance and Arrhenius parametrization 

We will focus on a situation when at each particular 
time moment the transition rates are in detailed balance, 
which means that they can be parameterized as 

fc,„ = fce-^(^=-^'), (4) 

where kia is the rate of a jump over the link a, starting 
with node i, and /3 = l/lksT), with ks and T being the 
Boltzmann constant and temperature, respectively, and 
k is the constant that insures a proper dimensionality. 
The choice of this constant does not influence the physics 
discussed below so, in what follows, we will set k = 1. 
The parameters Ei and Wa have the physical meaning 
of the potential well energies and the energy positions 
of the barriers, respectively [29, 30]. Thus, we will call 
parameter Ei the energy of node i and we will call Wa the 
barrier height at link a. The master operator depends 
on time H{t) = H{x{t)) via its dependence on the set 
X — {E, W) of parameters. 

In the detailed balance case the zero mode of the mas- 
ter operator H is usually referred to as the equilibrium 
distribution, and is represented by the Boltzmann distri- 
bution p^ , given by 

pf = Z-le-''^^ Z^^e-''^^ (5) 

k 

The non-zero eigenmodes ^ of H that describe relaxation 
(and are therefore referred to as relaxation modes) satisfy 
the zero total population condition = 0. Hence- 

forth, we will refer to the subspace of population vec- 
tors with zero total population as the space of relaxation 
modes or, simply, the relaxation space of H. Obviously, 
the relaxation space is an invariant subspace of H. 

We note that, strictly speaking, the expression for the 
rates should have a form kia = e~^(^°~^'^+'^'~'^° , where 
Si and Sa are the entropy terms corresponding to the 
metastable state, i, and transition state, a, respectively, 
however, for the sake of convenience the entropy terms 
will be absorbed by assuming that Ei is the free energy 
in state i etc. 

In order to simplify our notation, we will use parame- 
ters, Kj and ga, given by 



so that kja — ^jga' 

C. Driving protocols and generated currents 

Let Aix be the space of parameters x = {E, W), asso- 
ciated with a graph X. In this series of two manuscripts 
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FIG. 2. Link a between two sites i and j. According to our 
notation, da = {i,j}, doa — {i}, and dia — {j}. 



we study the currents generated by changing the system 
parameters with time so that at each particular time the 
system is in detailed balance. Thus, a driving protocol 
is a smooth map from time to the parameter space Aix- 
We will focus on periodic driving, which means that a 
driving protocol is a periodic function of time, with the 
values in Aix and time period td- Being specifically 
interested in the adiabatic limit we introduce the dimen- 
sionless time, r = t/T]j (here D is for "driving"), so that 
a driving protocol is fully characterized by its time period 
T]j and a periodic function s(r) with period 1. Obviously, 
s(t) can be conveniently interpreted as a map of a circle 
of unit length to Mx- The master equation adopts a 
form 



drP(r)=rz5H(s(r))p(r). 



(7) 



A stochastic (also known as empirical) current is rep- 
resented by a vector J = {Ja\ct £ A"i). For a stochastic 
trajectory rj the generated current is given by J = q/t, 
and Qa is the number of times the particle jumped over 
the link a in the direction of the orientation minus the 
number of times it jumped over the same link in the op- 
posite direction. 

The boundary map d can be extended to an operator 
that acts on vectors over the links into vectors over the 
nodes, with the physical meaning of being a discrete ana- 
logue of the divergence operator. By acting on the vector 
of currents, J, by our definition, it gives 



'9a={j,i) 



da={i,j) 



(8) 



aeXi 



where the summation, Yla'eXi'''^ ' should be understood 
as the sum over all links of the graph X having node i as 
the origin and pointing to the node j. 

We use the operator d because it allows us to write 
various expressions that involve current on a graph in an 
invariant form, i.e. in the form that does not depend on 
the choice of orientation of links of a graph. It is also very 
helpful, for the same purpose, to introduce an operator 
. The operator acts on distributions, resulting in 
currents. It has a physical meaning of a discrete analogue 
of the gradient operator (with minus sign), and is given 

by 

{d^p)a = Pdaa - Pdia for doa ^ dia, (9) 
and {d^ p)a = 0, when a is a self-link. 



It is straightforward to verify that the master operator 
can be recast in a form that resembles a Fokker-Planck 
operator for Langevin processes, which will turn out to 
be useful for our applications: 



dia.g{Wi, . . .W\Xx\) and E 



(10) 



where W = 

diag(i5i, . . . i?|Xo|) are diagonal matrices with val- 
ues of, respectively, barriers and node energies staying 
on the main diagonal. The stochastic current, generated 
over a time period [0, t^] by a driving protocol (not nec- 
essarily periodic), s(t), and averaged over the stochastic 
trajectories and initial distribution p(0), is given by 



J(t) = e-^^Wate^^«p(t), 



(11) 



where 



p{t) - Texp ^ dt'His{t'))j p(0) (12) 

is a formal solution of the master equation, and T is the 
time-ordering operator. We will focus on the long-time 
asymptotic behavior of the current generated via a peri- 
odic driving protocol. The average current, Q, defined 
as the average vector of fluxes through the graph links 
per period of the driving protocol, rather than per unit 
time is given by 



Q = N 



Ntd 



dtJ{t), N 



(13) 



In the N ^ oo limit, when substituting the general ex- 
pression for J [Eq. (11)] into Eq. (13) we can replace p{t) 
with the periodic solution of the master equation, which 
results in the following expression for the average current 
generated per period of the driving protocol 



Q(s;/3,td) = T£) / dr 
Jo 



(14) 



with p{t) being the periodic solution of the master equa- 
tion [Eq. (7)], written in terms of dimensionless time. 
In the present manuscript, we consider only the aver- 
age current Q(s;/3,T£)), generated over the driving pro- 
tocol period, with the focus on the low-temperature adi- 
abatic limit. Its fluctuations are studied in the second 
manuscript of the series. 



D. Good and bad parameters 

The quantization of currents, which we discuss below, 
occurs when the driving protocol belongs to a specific 
subspace of the parameter space. To define this subspace, 
we introduce the space Mx C Mx of good parameters 
for integer quantization. It is identified by the condition 
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FIG. 3. (a) Closed path along the contour 1 2 ■■•—>■ 6. 
(b-d) The collection of minimal spanning tree graphs repre- 
senting connected components obtained by successive elimi- 
nation of links with highest barriers. Orange color indicates 
the paths on the trees, lij{Xw), that contribute to the total 
integer- valued current. 

that for {E, W) £ Mx there is no simultaneous degen- 
eracy of the lowest node energies and a degeneracy of two 
barriers. It is convenient to describe Mx introducing 
a discriminant set, Dx, of bad parameters the protocol 
should avoid for the integer quantization effect to be in 
place, so that Mx = Mx — Dx- More formally, we say 
that {E, W) S Dx, if Wa — for some distinct edges 
a 7^ 7 (barrier degeneracy), and Ei = Ej < Ek for some 
distinct nodes i j and any other node k (energy degen- 
eracy). We can naturally represent Mx = UqUUi, with 
Uq and Ui standing for the sets of parameters {E, W) 
with no lowest-node energy, and no barrier height degen- 
eracy, respectively. 

Imagine now that we restrict our studies to protocols 
that keep some of the node energies permanently degen- 
erate and/or keep some of the link barriers permanently 
degenerate with each other. When studying such a set- 
ting with permanent degeneracy among some energies Ej 
and some barriers Wa it is convenient to introduce the 
sets Xo and Xi of the subsets J C Xo and 21 C of the 
sets of nodes and links, respectively, with permanently 
identical (degenerate) values of the energies and the bar- 
riers. Obviously Xq — Uj^x^Ji and X\ — UaeA-iSl. We 
denote Ej and W% the values of the (degenerate) ener- 
gies and barriers in J and 21, respectively. We can extend 
the definition of good parameters to this case and intro- 
duce the space Mx C Mx of good parameters for ratio- 
nal quantization, by forbidding simultaneous degeneracy 
among Ej and Wa- Similar to the integer quantization 
case, we have Mx ~ Uo^Ui. 



E. Spanning trees, associated integer currents, and 
Boltzmann distributions 

A tree is a connected graph without loops. A spanning 
tree of a graph is a tree subgraph that contains every 
vertex of a graph. Fig. 3 shows an example of a graph 
and three examples of spanning trees of this graph. 

For two distinct nodes of X and a spanning tree, X G 
X, we can identify the unique path, lij(X), that connects 
i to j through the spanning tree X. We also denote 



by Qij{X) the integer- valued current, associated with 
the described above path. More formally, a component 
Qij,a{X) — ±1, when the link a ~ {k,k'} belongs to 
the path lij{X) with the orienting arrow going k ^ k' 
and k' — )■ fc, respectively; and Qij,a{X) = 0, otherwise. 
For example, in Fig. 3(b), the orange color shows the 
path /i3(X(-b)) that generates a unit current through links 
(12) and (23). More generally, we denote by Q{1) the 
integer- valued current, associated with an oriented path 
I in our graph. Within the introduced notation we have, 
in particular Qij{X) = Q{lij{X)). 

We further introduce the energy W{X) — J^aeXi 
associated with a spanning tree X, and the Boltzmann 
distribution in the space of spanning trees 

4=Z'ie-^^W, Z = ^e-^^W. (15) 

X 

The integer- valued currents Qij(X) and the Boltz- 
mann distributions are important ingredients of the 
explicit expressions for the average currents, derived in 
this paper. 

F. Conserved currents, weighted scalar products, 
and non-circulating currents 

In this article we will frequently encounter conserved 
currents that are formally defined by the condition dJ = 

0, meaning that the sum of incoming minus the sum of 
outgoing currents at any node is zero. The pump current 
per cycle of the driving protocol (14) is an example of a 
conserved current. Indeed, if it were not true, then by 
repeating the cyclic protocol we would be able to increase 
the population of some nodes indefinitely. To avoid this 
inconsistency we have to assume that 

aQ(s;/3,Tz,)=0, (16) 

1. e. the sum of going through any node pump currents 
is zero. This condition restricts the number of possible 
independent pump currents to the number of indepen- 
dent loops of the graph. This number is also equal to 
the number of links that can be removed from the graph 
without breaking this graph into disjoint components. 

Another example of a class of conserved currents, im- 
portant for our applications, is represented by the cur- 
rents Q{1), associated with paths I in the case when I 
is closed, i.e., it is a loop. It is evident that Q{1) is a 
conserved current when I is a loop. We can say even 
more: the vector space of conserved currents is gener- 
ated by Q{1). A natural way to construct a basis set in 
the vector space of conserved currents starts with choos- 
ing a spanning tree X G X of our graphs. Each link 
a e Xi\Xi that does not belong to the tree gives rise to 
a loop laiX) obtained by concatenating the link a with 
the only non-self-intersecting path that connects its ends 
via the tree X. The set {la{X)\a G Xi \ Xi} forms a 
basis set in the space of conserved currents. 
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Conserved currents naturally do not re-distribute the 
population between the nodes. We will show that, when 
populations are re-distributed, the (non-conserved) cur- 
rents responsible for this process should minimize the 
circulating components (i.e., represented by conserved 
currents), which can be formulated as orthogonality to 
all conserved currents with respect to a weighted scalar 
product with conserved currents, defined by 

Hereafter, we refer to the currents orthogonal to the con- 
served counterparts, in the sense of weighted scalar prod- 
uct, given by Eq. (17) as non- circulating currents. Since, 
as stated before, conserved currents are generated by the 
currents, generated by the loops, the non-circulating (or- 
thogonality) condition for a current J reads 

(Q(Z), j)g ^ E QMgaT^Ja = 0, (18) 

aeXi 

for all closed paths (loops) in our graph. If we inter- 
pret ga and g^^ as the conductivity and resistance, re- 
spectively, of link a, followed by naturally interpreting 
{9a)~^Ja as the voltage on link a, then Eq. (18) claims 
that the sum of voltages of the links over closed path I 
is zero. Note that the factors Qa(l) = 0, ±1 restrict the 
summation in Eq. (18) to the links that belong to the 
path, as well as ensin-e that the voltages (potential dif- 
ferences) are always taken in the direction of the path, 
as they should be. 

G. Vector-potential and geometric properties of 
pump currents 

It has been demonstrated previously [29] that in the 
case when Arrhenius parametrization is fulfilled at all 
times, the pump current is totally determined by the ge- 
ometry of a contour in the space of population vectors, 
i.e., 

Q{s;l3,Tn)= [ Y.A'^dpk, (19) 

where the vector-potential. A, is a vector in the space 
of currents (i.e., the vector space spanned on the set Xi 
of the graph links) and a linear functional in the space 
of population vectors with zero total population. Stated 
equivalently, A is a linear operator that maps the popu- 
lation vectors with zero total population to the currents. 
We further note that, in the adiabatic limit, the formula 
(19) is almost obvious. Indeed, for a time-independent 
case the system is at equilibrium, represented by the 
Boltzmann distribution with zero average currents. 
When parameters are adiabatically time-dependent, the 
instantaneous average current, to the lowest order in 
time-derivatives, must be a linear combination of first 



time-derivative of the equilibrium node populations, i.e. 
J = A'^pf . Integrating over time of the protocol we 
find that in the adiabatic limit the total current is to- 
tally determined by the protocol geometry and does not 
depend on time explicitly, that is, 

Qis;p)= Vim Q{s;I3,td)= ly^A'^dp^. (20) 

■'^ k 

At this point we would like to note that, while A is de- 
scribed as a linear operator in an unambiguous way, the 
description of its components A'^ needs some clarifica- 
tion, due to restricted nature of the operator domain. De- 
noting by |j) the population vectors with pk = 6kj, and 
introducing a reference node a G Xo,we can see that the 
set of vectors |j) — |a) , with j G Xq and j ^ a form a basis 
set in the space of population vectors with zero total pop- 
ulation. We further naturally define A'^ = A{\k) — \a)). 
Obviously, the set (A.J|fc € Xq) of components depends 
parametrically on the choice of the reference node a, how- 
ever, A'^dpk is a-independent, since dpk = for 
normalized distributions. Therefore, the subscript a can 
be dropped in the expressions like Eqs. (19) and (20). 
Also the notation A'' — A^ is free of the above ambigu- 
ity, since — A'^ = A'^, are o-independent. 



III. MAIN RESULTS 

In this first manuscript of the series we study the 
average current Q{s; (3,td), generated during the pe- 
riod of a driving protocol in networks, represented by 
arbitrary graphs, with the main focus on its adia- 
batic followed by the low-temperature limit Q{s) = 
lim^_5.co linirD^oo Q(*; /5j "J"!))- Recently we have pro- 
posed the so-called Pumping- Quantization Theorem 
(PQT) [15] that describes the average currents generated 
per driving period, and can be formulated as follows: 

(i) Integer PQT: For an arbitrary graph and under 
the condition that the periodic driving protocol always 
stays inside the subspace of good parameters for integer 
quantization (which was defined in subsection II D), the 
average current per cycle must be quantized in the adi- 
abatic limit, which is followed by the low temperature 
limit. 

(ii) Fractional PQT: If we allow permanent degenera- 
cies of some barriers and/or energies and assume that the 
driving protocol belongs to the space of good parameters 
for rational quantization, while keeping all other condi- 
tions of (i) fulfilled, the current per cycle of the driving 
protocol is a rational (fractional) number. 

The main result of this first manuscript of our series 
is the averaging formula (AF) that quantifies the phe- 
nomenon of (fractional) quantization of currents. The 
AF states that the average current, Q{s), per cycle, s, 
of the driving protocol can be written explicitly in the 
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following form: 

^(«) = ;7- E (21) 

(fe,X) 

that expresses Q{s) as a normalized sum of the integer 
valued conserved currents, Q{s;k,X), associated with 
(fe, X), which represents all possible ways to split spe- 
cific parameter degeneracies that are encountered along 
the driving protocol, referred to as global degeneracy res- 
olutions, with the normalized integer denominator nx 
that represents the number of legitimate global degener- 
acy resolutions. The explicit construction for Q{s; k, X) 
and definitions of nx and {k,X) are provided in sec- 
tion VII. At this point we would like to emphasize that 
the AF for the average current not only delivers a proof 
of the fractional PQT, but also provides an explicit con- 
struction that allows the average current to be computed 
in the low-temperature adiabatic regime. Note that in- 
teger quantization can be viewed as a particular case; of 
rational quantization when there is no degeneracy, which 
implies nx = 1. We will first consider in section VI 
this special case of integer quantization in some detail, 
while the more general case of rational quantization is 
postponed to section VII. There are two main reasons 
for that. First, the integer- valued currents Q{s;k,X) in 
the r.h.s. of the AF [Eq. (21)] can be viewed as integer- 
quantized currents that correspond to legitimate ways, 
labeled by {k,X), to resolve the degeneracy, so that un- 
derstanding the nature of integer quantization provides 
an additional insight into the rational quantization phe- 
nomenon. Second, at the present point the nature of 
integer quantization is understood on a deeper level of 
the current probability distribution, or more specifically 
the associated with it large deviation function S{Q), as 
presented with some detail in the second manuscript of 
the series. 

On the way to the derivation of the AF, we will obtain 
several results that are of independent interest in the the- 
ory of the pump effect. Thus the derivation of the AF 
[Eq. (21)] includes the following steps: 

1) Wc show that components of the vector-potential, 
A in Eq. (19) can be identified with the current distri- 
butions on a graph of an electrical network, in which the 
links represent the resistances r^, given by 

r-a=5a'=e^'^% (22) 

so that the component, A'^, corresponds to current vector 
in our network that would appear if a unit current enters 
the network and leaves it at nodes a and k, respectively. 

2) Wc show [Eqs. (25) and (26)] that the identifica- 
tion of A with a current in a resistance network allows 
us to apply the Kirchhoff's theorem and determine the 
expressions for A'^ and A'' — A'^ explicitly in terms of 
the Boltzmann weights of all spanning trees and currents 
that correspond to paths that connect nodes k and k' 
on these trees. More specifically, we introduce an energy 
function on the set of spanning trees that associates with 



a spanning tree X the energy W{X) given by the sum of 
the barrier heights of the links that belong to X, as well 
as the corresponding Boltzmann distribution , given 
by Eq. (15). 

3) The above identification leads to a general and 
explicit expression [Eq. (27)] for the average current 
Q{s; I3,T£)) in terms of the tensor Q*^"^ (with values in 
the integer-valued currents, and superscripts k and X 
running over the graph nodes and spanning trees, respec- 
tively) that depends on the graph structure only, and a 
distribution 7 in the combined space of nodes and span- 
ning trees. Explicit integral expressions [Eqs. (28) and 
(30)] for 7 involve only two loops (closed paths) in the 
spaces of node and spanning tree distributions: these are 
the periodic solution /?(t) of the master equation and the 
time dependent tree Boltzmann distribution Q^{t) that 
depends on time via periodic dependence W(t) of the 
barrier data, according to the driving protocol. 

4) The adiabatic td — 00 limit is obtained in a 
straightforward way by replacing p{t) in the expressions 
for 7 [Eqs. (28) and (30)] with its adiabatic limit, rep- 
resented by p^(t). The obtained expressions [Eqs. (32)] 
show symmetry between the sets of nodes and spanning 
trees, which might be indicative of some kind of duality in 
the adiabatic limit, the latter being not yet understood. 

After establishing the AF, we develop a winding in- 
dex interpretation of the average current. This allows 
the space Mx of robust parameters, or simply the ro- 
bust space to be identified, which constitutes the broad- 
est subspace that provides robust behavior of the pumped 
currents. This means that, for a periodic driving proto- 
col s that stays at all times within the subspace of ro- 
bust parameters, the limit lim^_5.oo IvxIt-^^oc Q{s) does 
exist, and the generated current does not change upon 
any deformation of the driving protocol s as long as s 
stays within the robust subspace. In particular, this 
means that the generated current Q{s) depends only 
on the topological class of the driving protocol, con- 
sidered as a 1-dimensional cycle in the space of robust 
parameters. It is convenient to represent the subspaces 
Aix = -Mx — Dx as a result of withdrawing the dis- 
criminant subset, Dx, of non-robust parameters. The 
discriminant set has codimcnsion 2, so it is possible to 
describe the current in terms of how many times the 
driving protocol winds around the components of the dis- 
criminant set Dx- We explore the structure of Dx and 
show that it consists of a bunch of (codimension 2) cells, 
glued together via their borders, and is closely related 
to the structure of the discriminant Dx of bad param- 
eters that was introduced in subsection II D. The only 
difference is that certain simultaneous degeneracies do 
not contribute to the currents and, therefore the corre- 
sponding cells can be withdrawn from the discriminant 
set, resulting in Dx C Dx- 

Since the discriminant subspace Dx, and hence Dx, 
has codimension 2, robust driving protocols, i.e., the ones 
that stay within the robust parameter subspace, are typ- 
ical. This means that, if one creates a protocol with 
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randomly chosen parameters, such protocol will belong 
to the robust space. 



IV. EXPLICIT EXPRESSIONS FOR AVERAGE 
CURRENTS AND KIRCHHOFF THEOREM 

A. Vector potential and Kirchlioff problem 

In this subsection we derive the expression for the av- 
erage pumped current, given by Eq. (19) by (a) express- 
ing the current J, defined by Eq. (11) as a solution of 
a linear problem, (b) demonstrate that the latter has a 
unique solution, which allows us to represent the current 
in a form J = A'^pk, (c) We will further demonstrate 
the equivalence between the aforementioned problem and 
the Kirchhoff problem of finding the current distribution 
in an electrical network with the given distribution of the 
(incoming and outgoing) external currents. This analogy 
will allow us to apply the Kirchhoff theorem that provides 
an explicit form for the sohition of the Kirchhoff prob- 
lem, and represents the linear operator A as a weighted 
sum over all spanning trees of X. This will be presented 
in subsection IV B. 

The linear problem, discussed above is formulated in 
terms of two properties of the generated current J: 

(i) the continuity equation dJ = p, 

(ii) J is non-circulating, i.e., it is orthogonal to all 
conserved currents with respect to the weighted scalar 
product, i.e., (J, Jo)g = for any Jq with dJo = 0. 

To establish property (i) we apply the divergence oper- 
ator to the definition of J, given by Eq. (11), followed by 
making use of the representation of the master operator, 
given by Eq. (10), which results in dJ = Hp = p. Prop- 
erty (ii) is verified by the following sequence of equali- 
ties based on the expression for J [Eq. (11)] and defi- 
nition of a conjugate operator (J, Jo)s = e'^^J • Jq = 
d^e^^p ■ Jo = el^^p ■ dJo = 0. 

Provided the linear problem under consideration has a 
unique solution, it can be represented in a form J = Ap, 
where j4 is a linear operator that acts from the vector 
space of relaxation modes to the vector space of cur- 
rents. This immediately establishes the expression for 
Q{s), given by Eq. (19), as well as fully describes the 
linear operator A as providing the solution of the linear 
problem under discussion. 

The uniqueness of the solution of our linear problem 
is almost obvious. Let J and J' be its solutions. Then 
{J— J') is a conserved current orthogonal to all conserved 
currents, and therefore is equal to zero, i.e., J = J' . The 
existence of a solution follows from the facts that the 
boundary (discrete divergence) operator d maps the vec- 
tor space of currents onto the relaxation space of H, and 
the kernel of d coincides with the subspace of conserved 
currents. Therefore, since the space of non-circulating 
currents is the orthogonal complement to the space of 
conserved currents, the operator d provides an isomor- 



phism between the space of of non-circulating currents 
and relaxed populations. Therefore A can be chosen as 
the inverse of the above isomorphism. 

To establish equivalence between our linear problem 
and the Kirchhoff problem, we recast the conditions (i) 
and (ii) in terms of the components A'^, i.e., setting pj = 
Sj — Sj, so that condition (i) reads 

{dA',)j=Sl-5]. (23) 

We further recast condition (ii), i.e., the orthogonality 
condition using conserved currents, associated with loops 
I of our graph [Eq. (18)] 

J2 Qa{l)ef'''--Al^=Q. (24) 

We now interpret as electrical current that goes 
over link a of an electrical circuit, represented by our 
graph, provided we have an incoming and outgoing cur- 
rents both of value -|-1 at nodes a and k, respectively. 
According to the discussion that followed Eq. (18), we 
interpret = e^^^" as the resistance of link a, which 
makes e^^°A'l^ to be the voltage on link a. Within the 
presented interpretation, Eq. (23) represents first Kirch- 
hoff law: At any node in an electrical circuit, the sum 
of cmTcnts flowing into that node is equal to the sum of 
currents flowing out of that node; whereas Eq. (24) re- 
flects second Kirchhoff law: The directed sum of the elec- 
trical potential differences (voltages) around any closed 
circuit is zero. Therefore, is given by the solution of 
the Kirchhoff problem: Find the stationary currents that 
flow over the links of an electrical network, represented 
by a graph X, provided we have an incoming and outgo- 
ing external currents both of value -|-1 at nodes a and k, 
respectively. 

B. Kirchhoff theorem and explicit expressions for 
average currents 

The solution of the Kirchhoff problem is provided by 
the celebrated Kirchhoff theorem [33, 34]. The Kirchhoff 

theorem provides an explicit expression for the current, 
in our case denoted A^, in a form of a weighted sum 
of integer currents Qak{^) (defined in subsection HE) 
over all possible spanning trees X G X oi our graph 
[33, 34]. According to Kirchhoff 's theorem, the relative 
(non-normalized) weight of a spanning tree X is given by 
the product of conductances (inverse resistances) of its 
links, the normalization factor ensures the (normalized) 
weights to sum up to one. Given the above interpretation 
of e^^° as the resistance of link a, the definition of the 
tree energy W{X), and the tree Boltzmann distribution 
[Eq. (15)] Kirchhoff's theorem leads to an explicit 
expression for the vector potential A 

Al = Y,elQak{x), (25) 

X 
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which also implies 



(26) 



X 



The explicit expression [Eq. (26)] is very convenient for 
obtaining the low-temperature limit /3 — >• oo, which will 
be explored in the manuscript. Wc wish to reiterate that 
in Eqs. (25) and (26) summation runs over all spanning 
trees X (Z X of our graph X. 

An explicit expression for the average current is ob- 
tained upon substitution of Eq. (25) into Eq. (19). To 
emphasize some symmetry between the nodes and span- 
ning trees we also introduce a notation = Qak(X), 
with omitting a, when the expressions are independent 
of the reference node choice. This results in 



k,X 



with 



%xis;P,rD) 



Qxdpk- 



(27) 



(28) 



We can think about 7 as a distribution in the combined 
space of nodes and spanning trees. The distribution 7 
satisfies the following properties: 



(29) 



The first property allows the reference node to be omitted 
in the r.h.s. of Eq. (27), the second property provides the 
conserved nature of the average current Q{s). 

Integration by parts and the Stokes theorem provide 
two alternative expressions for 7 



7fcjf(s;/3,'rD) 



PkdQ% = / dg%Adpk, (30) 



where A is surface spanned on s in the cartesian product 
of the space of normalized node distributions and the 

space of normalized tree distributions. Since both dis- 
tribution spaces are contractible, such a surface always 
exists. 

The derived expressions for the average current are 
very convenient for deriving its adiabatic limit Q{s;t£)), 
which is achieved by simply replacing -^[s] I3,td) in 
Eq. (27) by its adiabatic limit 7(5, /3), which boils down 
to replacing the distributions p in the expressions, given 
by Eqs. (28) and (30) with the Boltzmann distributions 
p^, which results in 



g(s;^) = ^g'=^7,x(«;^)> 



(31) 



with 



Ikx 



k,X 



B J„B 



Qxdp, 



= / dg%Adp^. (32) 



We note that the adiabatic expression for the distribu- 
tion 7 show symmetry between the sets of nodes and 
spanning trees of the graph. This might bo a signature 
of certain duality, with the nodes of the dual graph being 
represented by the spanning trees of the original graph 
X. At this point it is unclear to us whether such duality 
really takes place, and in particular how to define the 
links of the dual graph that connect the spanning trees. 



V. TOY MODEL 

As a simple example of a graph of the type shown 
in Fig. 1, we consider a model illustrated in Fig. 4: a 
particle can randomly jump betwecni two sitcis along two 
different paths. Each path can be traversed in both di- 
rections and arrows on the links indicate the directions of 
particle transitions that contribute to the currents with 
the plus sign. This is the minimal model that shows in- 
teger current quantization and that we will work out in 
this section in detail. A reader, who is interested to see 
further specific examples that include fractional quanti- 
zation, might want to look at our previous publication 
[15]. 



A. Markov-chain on a 2-nodes-2-links graph 

Populations on the nodes in our toy model are de- 
scribed by the population (probability) vector, p = 
iPi,P2), where pi and p2 are probabilities for a parti- 
cle to be in the first and the second nodes, respectively. 
They evolve according to the master equation [a partic- 
ular case of Eq. (1)], 



d^ 
dt 



H 



(33) 



with H, being the master operator, according to its gen- 
eral definition [Eq. (2)]. In terms of the parameters ga, 
with a = 1, 2, and kj, with j = 1,2, the master operator 
H reads 



H 



-Ki(.9i 



^.92) 

92) 



1^2(91 ^ 
-K2{gi 



92) 
^92) 



(34) 



B. Driving protocols and generated currents 

For the model in Fig. 4, the 4-dimensional parame- 
ter space M describes the sets x = {Ei,E2,Wi,W2) 
of controlled parameters. A driving protocol s is nat- 
urally a cycle in M. For example, s can be a cy- 
cle defined by E2 = W2 = 0, E-i{t) = Eocos{2nT), 
Wi{t) = Wq sin(27rT), where t is the dimensionless time 
given by r = I/td', with td and t being the period of 
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FIG. 4. The "2-nodes-2-links" model of stochastic transitions 
between two states along two possible paths. Different paths 
are characterized by the values of the corresponding potential 
barriers Wi and W2 ■ The two states are characterized by the 
values of their well depths, Ei and E2. Periodic modulation 
of these parameters leads to a nonzero, on average, circulat- 
ing current in the preferred (clockwise or counterclockwise) 
direction. 

driving and actual time, respectively. The driving proto- 
col is, therefore, completely determined by its duration 
td and s{t) for r G [0, 1). 

In addition to node populations, pi and p2, we consider 
the vector of currents J = (Ji, J2), with Ji and J2 being 
the random variables that describe how many times per 
unit time the particle jumps through the corresponding 
graph links, as shown in Fig. 4, counting counterclockwise 
and clockwise transitions with a plus and minus sign, 
respectively. The boundary operator d, which returns 
the sum of currents outgoing from each node when it is 
applied to the vector of currents, takes the form 

'-{:\-\)^ <-> 

in the full agreement with its general case definition 
[Eq. (8)].^ 

Its conjugate operator returns the difference of node 
populations between the origin and the target nodes of 
each link, namely. 

By introducing operators, W = diag(Wi,M^2) and E = 
dia.g{Ei, E2), so that 

we can express the master operator, in the form, given 
by Eq. (10). 

C. Geometric properties of vector-potential A in 
toy model 

The vector potential A can be derived in adiabatic 
limit by considering small deviation of the probabil- 
ity distribution from its instantaneous equilibrium, p — 



-f 5p, where p^ = Z-^el^^' ,e^^^), Z = {e^^' + 
gPE2y Substituting this into the master equation (33) 
and disregarding higher time-derivatives than the first 
one, we can find 6p as a function of dp^ /dt, i.e. 

-dpf/dt 

dpi = -dp2 = 7 ■ r-, ■ r. (38) 

(51 +52)('«1 + H2j 

We then find the average current using Eq. (11), which 
leads us to 

•h = giKiSpi - giK26p2, 

(39) 

J2 = -92HlSpi +g2K2Sp2- 

Due to Eqs. (38) and (39), the differential of the total 
current, dQ — Jdt, can be written as 

dQ = A^dp^+A^dp^, (40) 

where A^ = a{gi + 52)"^(-ffi, 92), = (a - 1)^^, 
and where a is arbitrary constant, which has no physi- 
cal meaning. Its arbitrariness reflects the fact that, due 
to particle conservation, the vectors A-^ and A^ are not 
independent. It is rather their difference, 

A^-A' = {g,+g2)-'igu~g2) (41) 

that can be uniquely defined. 

To make parallels with the more general Eq. (26) , note 
that we can identify A^ — A^ by introducing the vector 

= (5i+52)~'(.9i, .92) (42) 

of Boltzmann weights with components on both of the 
graph's spanning trees, shown in Fig. 5. We can further 
introduce the vectors of currents, 

QuiXi) = (1, 0), Qi2{X2) = (0,-1), (43) 

which are generated by the paths that start at site 1 
and end at site 2 staying, respectively, on the tree Xi in 
Fig. 5(a) and on the tree X2 in Fig- 5(b). Then, 

A^-A'= J2 QM)q%. (44) 

Xi;i=l,2 

which shows that Eq. (41) is a special case of Eq. (26). 

We would like also to discuss the analogy between 
the properties of the vector potential A and Kirchhoff's 
laws for electric circuits. The components of the vector 
A^ — A} can be viewed as currents on the two links of 
the graph. An important observation is that those cur- 
rents are the same as they would be if we consider the 
currents in a very simple electric problem, in which the 
links represent the resistances R\ — gi^ and R2 — g2^ 
for the first and the second link, respectively, and the 
nodes represent the junctions where a unit current en- 
ters (first node) and leaves (second node) the network, 
respectively. According to the Kirchhoff's laws, the drop 
of the voltages around the loop of the network must be 
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w. 



(b) 



El 



FIG. 5. Minimal spanning trees for (a) Wi < W2 and (b) 
Wi > Wi. At the low temperature limit, motion is restricted 
to be on a tree graph that, for toy model, consists only of 
one link and two nodes. The choice of the minimal spanning 
tree at each time interval depends on the relative strengths 
of barriers, e.g. interval (iv) corresponds to the tree (a) and 
interval (ii) corresponds to the tree (b). 



J,„=l 




FIG. 6. Resistance network that corresponds to the model in 
Fig. 4. 



zero, whereas the sum of the currents entering any node 
must be zero, which can be expressed as 



Jii?i + J2R2 = 0, 
1 + J2 - Ji - 0. 



(45) 



The solution of the above linear system [Eq. (45)] yields 
(Ji, J2) = (51 + 52)~"'^(ffij —92), which coincides with 
Eq. (41). So we can identify the components of the vector 
J that solves the Kirchhoff's laws with components of 
the vector-potential, — , e.g., the latter also solves 
Eq. (45). Note also that by using our notation, after we 
replace (Ji, J2) by the components of — A^, Eq. (45) 
can be compactly rewritten as 

{A^~A\Jo)g = 0, diA' - A% ^ Sf - dl (46) 

where Jo = (Ij 1) is a unit conserved current that circu- 
lates around the loop in Fig. 4. 



D. Current quantization in the low-temperature 
limit 

For our toy model, the space Dx of parameters with 
bad degeneracies is defined, simply, by the conditions 
El = E2 and Wi = W2- Any driving protocol that avoids 
this degeneracy space, hence, belongs to the space M.x 
of good parameters. For example, consider a protocol, 
given by E2 = W2 ^ 0, Ex(t) = cos(27rT), Wx{t) = 
sin(27rr). One can distinguish four time segments with 
distinct dynamical behavior patterns: 



(i) for T G where e < tt/2, the path goes 
through the point Wi = W2 but we always have Ei > E2; 

(ii) for T e [e, TT — e], we have Wi > W2; 

(iii) for T G [tt — e^TT + e], we pass through the point 
Wi — W2 but we always have Ei < E2; 

(iv) for re [n + 6,271 — e], we have 14^2 > Wi. 

Consider stochastic dynamics on each of those seg- 
ments in the low temperature limit. Whenever the bar- 
riers, Wi and W2 are different, when p ^ 00, we can 
omit in Eq. (41) the terms e'^^^^ or e'^'^^ if Wi > W2, 
or W2 > Wi, respectively. This means that at the seg- 
ments (ii) and (iv) A^ — A} is given by (0, —1) and (1, 0), 
respectively. At segments (i) and (iii) we encounter a sit- 
uation when W\ = W2. However, at these segments we 
have Ex > E2 and E2 > Ei, which in the /3 — 00 limit 
means that the whole population is almost completely 
accumulated in the node 2 and 1, respectively. There- 
fore, the population vector practically does not change 
at these intervals, and, therefore, no significant average 
current is generated, according to Eq. (40). Stated dif- 
ferently the current is generated at the segments (ii) and 
(iv) only. Note also that at the segments (ii) and (iv) the 
population difference pi — p2 changes from -1 to 1, and 
from 1 to -1, respectively. 

Substituting the obtained values of A for the segments 
(ii) and (iv) into the integral expression for the average 
current [Eq. (20)] wc find that current generated by the 
protocol (i)-(iv) is given by 



Q{s) - (1,1), 



(47) 



i.e. the particle passes overall exactly one time through 
both links of the graph in the counterclockwise direc- 
tion. Note also that this current is a conserved current 
in the sense that sum of its components entering any node 
equals the sum of components leaving this node. 

The low- temperature average current, given by 
Eq. (47), can be interpreted using the following deter- 
ministic picture. If the barrier through the second link 
is higher, this link cannot be traversed and the motion is 
effectively restricted to the minimal spanning tree graph, 
as shown in Fig. 5(a). Similarly, when the first barrier is 
higher, the particle motion is restricted to the tree graph 
shown in Fig. 5(b). 

At the start of the protocol (i)-(iv), the particle distri- 
bution is concentrated at site 2 and all currents are ex- 
ponentially suppressed. When energy Ei becomes lower, 
the particle moves to site 1 via link 2, and still stays 
in site 1 while barrier 1 becomes preferable for a tran- 
sition and, therefore the particle returns to site 2 via 
link 1, resulting in quantized current given by Eq. (47). 
This picture, of course, is oversimplified because we apply 
the adiabatic limit before the low-temperature counter- 
part, so that the true motion is stochastic rather than 
deterministic. In fact, a particle can make, and actually 
makes a large number random transitions through both 
links during the cycle duration and quantization that we 
discuss appears only on average. Nevertheless, the deter- 
ministic picture correctly predicts the average current. 
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E. Winding index 

The parameter space Mx for our toy model is 4- 
dimensional. The space Dx of bad parameters is 2- 
dimensional, being defined by two constraints, Ei = E2 
and W2 = W2- To visualize the structure of Dx^ we 
fix one of the barrier heights, W2 = 0. The space Dx of 
bad parameters is represented by a straight purple line in 
Fig. 7. In the adiabatic low-temperature limit any driv- 
ing protocol that encloses this line produces an integer- 
valued current, given by the protocol winding number 
around this line, i.e., the number of times the curve 
in the good parameter space that represents the proto- 
col winds around the straight line that represents Dx- 
There is an analogy between the above topological rep- 
resentation of pump currents and the Aharonov-Bohm 
effect in quantum mechanics. In the latter, the phase 
of the electronic wave function changes upon enclosing a 
quasi-one-dimensional solenoid with a magnetic field by 
an amount proportional to the total flux of the field in- 
side the solenoid. This analogy follows from the similar- 
ity of the expression for pump current and the geometric 
(Berry) phase in quantum mechanics that has been dis- 
cussed previously in the literature [1, 11, 14, 35, 36]. 

To make the above analogy more explicit, as well as 
to understand better the nature of the low-temperature 
limit we consider a contour that encloses the bad param- 
eter space only once, as shown in Fig. 7. By applying 
the Stokes theorem to Eq. (20), we can express the cor- 
responding pump current, Q{s), generated at arbitrary 
finite temperature, in terms of an integral over a surface, 
A, spanned on the contour, s that represents the driving 
protocol 



Q{s) = / dAAdp (48) 
J A 

= ^ j F°'\E,W)dEj AdW^, (49) 



where, in our toy model, z, a = 1,2, and we define the 
curvature F of the vector potential A in a standard way 



fc=l,2 



dE, 



(50) 



Using the mentioned above analogy with the phase of a 
quantum wavepacket we can view the curvature as a 
"magnetic field" in the parameter space [15]. For low 
enough temperatures this magnetic field is concentrated 
along the bad parameter space, as shown in Fig. 7. The 
total flux of this fleld equals 1. Decreasing temperature 
in the model does not change the total flux, while forcing 
the field to be concentrated in a narrow tubular neigh- 
borhood of the straight line that represents the bad pa- 
rameter space. In the /3 — >■ 00 limit this field is restricted 
to an infinitely narrow ID solenoid along the bad param- 
eter space. This leads to the above mentioned analogy 
with the Aharonov-Bohm effect. 




FIG. 7. Space of the parameters with bad degeneracies (pur- 
pur) and a contour in the space of good parameters (blue) 
that winds around the component of the bad parameter space. 
The winding number that corresponds to this contour is 1. 
Parameter W2 is set to zero. 



VI. INTEGER QUANTIZATION 

In this section, we return to the arbitrary graph ge- 
ometry and study the low-temperature adiabatic limit 
Q{s) — lini^_j.oo limru^oo Q(s; td) in the situation 
when the driving protocol stays within the subspace Aix 
of good parameters for integer quantization (see subsec- 
tion II D for the definition of A^x)- Let us introduce no- 
tation, k® X, which means the unit vector in combined 
population-current space. This vector has unit compo- 
nents on the node k and links of the tree X and has 
zero other components. The sum-over-trees expression 
for the adiabatic limit Q(s; /3) of the average current 
[Eqs. (31) and (32)] is very convenient for taking the 
low-temperature limit, which boils down to identifying 
the low-temperature limit of 7 = Ikx^ ® ^ 



Q{s)=J2Q'''%x{s),%xi^ 

k,X 



1™ 7fejf(s;^)(51) 

p— >-oo 



provided the /? — J' 00 limit in Eq. (51) does exist. In this 
section we explicitly calculate the low-temperature limit 
7(5) in the case when the driving protocol fully lies in the 
space A^x of good parameters for integer quantization, 
and show, in particular, that the components 7^,^^ are 
all integer, which implies integer quantization of average 
current in the setting under consideration. 



A. Maximal-weight and minimal spanning trees 

In the adiabatic limit, the distribution 7(5) is explic- 
itly expressed in terms of the Boltzmann distributions 
and in the node and spanning tree spaces, respec- 
tively. In the case when there is no degeneracy of the 
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lowest node energy Ek the Boltzmann distribution has 
a well-defined low-temperature limit = Sjk, i.e. the 
population is concentrated in the lowest-energy node k. 
The low-temperature limit of the tree Boltzmann distri- 
bution under the condition of no degeneracy between 
any two barrier heights is less obvious. The limit can still 
be identified by introducing a notion of the minimal span- 
ning tree Xw, and its two properties: (i) the tree Xw 
is non-degenerate, and minimizes the energy W{X), (ii) 
the tree Xw is totally determined by the ordering of the 
barrier heights, the actual values of the barrier heights 
within a given ordering is irrelevant. Property (i) com- 
pletely identifies the low temperature limit = Sj^-^^ 
of the tree Boltzmann distribution in the case of non- 
degenerate barrier data. Property (ii) justifies the nota- 
tion where [W] denotes the equivalence class of 
the barrier data with the same ordering of the barrier 
heights. 

The minimal spanning tree Xw of X, associated 
with non-degenerate barrier data W is a subgraph, con- 
structed by sequentially removing the links with the high- 
est possible barriers such that the remaining graph re- 
mains connected. Eventually, this process ends up with 
a tree (i.e. a graph without loops), referred to as the 
minimal spanning tree. 

The minimal spanning tree has an obvious, yet very 
important property that can be actually used as its al- 
ternative definition, and, therefore, referred to as the 
characterizing property of a minimal spanning tree: For 
each withdrawn link (i.e., a link that does not belong to 
Xw) Cf with da = {i,j}, the barrier Wa is higher than 
any of the barriers associated with the links of the path 
lij{Xw)- In particular, the physical meaning of the in- 
troduced current vector, Q{Xw) that corresponds to a 
path path lij{Xw)^ is that it is generated when a particle 
moves from node i to node j on the graph X, avoiding 
the highest barriers. 

Property (ii) immediately follows from the above ex- 
plicit construction of Xw that uses the ordering of the 
barrier heights only. Property (i) is demonstrated in ap- 
pendix A. In this subsection we just present a physical 
argument in support of property (i). According to the 
sum-over-trees representation the population redistribu- 
tion in the /3 — > oo is restricted to the spanning tree of 
the maximal weight, i.e., of the minimal value of W{X). 
On the other hand, random motion of the particle tends 
to avoid the links with the highest barriers, and, there- 
fore, should be restricted to the minimal spanning tree 
Xw, which supports the statement of property (i). 



B. Good parameters, minimal spanning trees, and 
quantized currents 

Current generation in the low-temperature adiabatic 

limit can be interpreted using a simple and clear in- 
tuitive picture. When the protocol goes over the sub- 



space Uq C A4x, when the lowest node energy is non- 
degenerate, the particle sits in the uniqiie lowest-energy 
node, and no current is generated at all. Current can be 
generated only during time intervals when at least two 
lowest node energies are close to their degeneracy. How- 
ever, in this case, by the definition of good parameters, 
the barrier data W is non-degenerate (i.e., the protocol 
goes through Ui C Mx), and population redistribution 
between any two nodes i and ), to avoid the highest barri- 
ers, is restricted to the tree Xw, and, therefore, produces 
the current proportional to Qij{Xw)- 

To make use of the above argument, we partition the 
time period of a driving protocol into a set of small 
enough segments, so that for each segment we either 
encounter only the minimal well depth degeneracy or 
just the potential barrier degeneracy (this can always be 
achieved due to the well-known argument of van Kam- 
pen). We refer to them as to 0- and 1-segments, respec- 
tively. If necessary, we merge the consecutive segments 
of the same type to make the 0- and 1-segmcnts alternat- 
ing. In the /3 — > oo limit, no current is generated at any 
1-segment, since the populations is concentrated at the 
node with the lowest energy value, while nonzero average 
pump current can occur only when the state probability 
vector changes [29]. Since the 1-segment ends also be- 
long to 0-segments, in the /3 — ^ oo limit the populations 
at the beginning and at the end of a 1-segment, labeled 
a, are concentrated in well-defined nodes, denoted by 
and k'^, respectively, and determined by the neighboring 
0-segments. 

As already argued above, when the particle moves from 
one node to another to generate current, it avoids the 
highest barriers i.e. it moves along Xw, which is asso- 
ciated with non-degenerate barrier data W . Since the 
paths with non-vanishing probabilities from fc^ to k'^ be- 
long to a tree Xw^, the total current Q", passing during 
the time of the 1-segment a has values 1 on any link 
on the shortest path Ik^k' {Xwa) and zero values on the 
other links. Finally, the current per cycle is gcncirated by 
the concatenation of the consecutive paths h^k'^iXw^) 
that correspond to the a-segments. This explicitly iden- 
tifies the generated integer-valued average current per 
cycle, 

Q{s) = ^Qk^k'MwJ, (52) 

a 

where the summation runs over all 1-segments, and, com- 
bined with the definition of Qij{Xw), Eq. (52) yields 
integer-valued Q{s). Since, by construction, the 1- 
scgments are alternating with the 0-segments, we have 
for consecutive 1-segments A;^ = ka+i, which implies 
dJ^a.Qknk'^i-^Wa) = 0' or equivalently Q{s) is a con- 
serving integer- valued (quantized) current. 

We have derived the expression for Q{s) [Eq. 52)] us- 
ing plausible and transparent physical argument. We are 
now in a position to sketch a rigorous proof of the integer 
PQT, or more specifically of the validity of Eq. (52). The 
sum-over-trees formula [Eq. (51)] allows us to identify 
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the low-tempcraturc hmit, 7(s), of the joint distribution 
7. We start with a formal notation for the aforemen- 
tioned partitioning of the driving protocol time period 
[0,1]/{0,1} = 51 = UaUa where /„ = [<_i,r„] 

and la = [Ta,T'g\ are the alternating 0- and 1-segments, 
respectively. We further recast the expression for the adi- 
abatic, still finite-temperature limit of the joint distribu- 
tion [Eq. (32)] as a sum of the segment contributions 

7(s;/3) = ^ / ° (ire^®p^ 



Xw = Property (iii) is established by represent- 

ing 



W ■ dwQ 



(56) 



(53) 



and apply integration by parts to the second term (sum 
over 1-segments), resulting in 

a 

-^e''(s(ra);/3)®p^(5(Ta);/3) 

a 

+ rfre^(s(r);/3)®p^(s(T);/3) 



[where W is considered as a vector in the space over 

the spanning trees with the components Wj^ = Ty(X)], 
and noting that E, W, p^ , and are bound, and the 
gradient of a Boltzmann distribution tends to zero in the 
/3 — > 00 limit in the case of non-degenerate lowest-energy 
state (node or tree) . The last property can be verified by 
obvious and straightforward estimates. 

We conclude this subsection with emphasizing that the 
presented derivation not only demonstrates integer quan- 
tization of pumped currents in the low-temperature adi- 
abatic limit, but rather provides an algorithmic way of 
computing the integer-valued current, given the driving 
protocol. The integer quantization formula (52) is, obvi- 
ously, a special case of the AF (21). 



/ dTQ^{s{T); igi p^{s{t);P). (54) C. Deterministic picture for integer quantization 



Before applying the low-temperature limit to Eq. (54) 
we would like to note that the expression is exact, and 
is valid for any partition of the time period into a set 
of alternating segments. In particular, it will be used 
in subsection VII C to identify the low-temperature limit 
for the rational quantization setting (permanent degen- 
eracies) . 

The /3 — )■ 00 limit of Eq. (54) for the integer quanti- 
zation setting under consideration is obtained using the 
following facts: (i) The parameter sets that correspond to 
the segment ends s(Ta),s(T^) G C/q H [/i represent non- 
degenerate lowest-energy and barrier data, and, there- 
fore, the node and tree Boltzmann distributions in the 
contact (the first two) terms in Eq. (54) have well-defined 
limits with the populations concentrated in the lowest- 
energy node and minimal spanning tree, respectively; (ii) 
the minimal spanning tree XvK(r) is the same for all times 
that belong to a 1-segment /«, so that a notation Xg. 
makes perfect sense; (iii) the integral terms (the last two) 
in Eq. (54) vanish in the (3 ^ 00 limit. Based on these 
properties the low-temperature limit of the joint distri- 
bution is identified as 



(55) 



and Eq. (52) is obtained by applying the bilinear form 
Q with the components Q'^'^ to Eq. (55), making use of 
the property Q{X (g) {k' - k)) = Qkk'{X). 

Property (i) is obvious. Property (ii) follows from the 
connected nature of a segment topology and the fact that 



In this subsection we illustrate the construction of inte- 
ger quantized current Q{s) as a sum of piecewise integer 
currents Qk^k' {Xwa)^ defined on minimal spanning trees 
(see Fig. 3). The path in the control parameter space can 
be split into overlapping intervals. Intervals of the first 
type have an unchanging node with the lowest energy, 
while intervals of the second type correspond to transi- 
tions between nodes but having the minimal spanning 
tree unchanged. 

For a graph in Fig. 3(a), we illustrate the determinis- 
tic evolution for some periodic protocol. Orange arrows 
indicate the path that dominates the average currents at 
low temperature. A particle starts at the minimal node 
1 having the largest barriers along the links connecting 
nodes {3,4}, {1,6} and {2,4}. By the time when this 
configuration of largest barriers changes the node with 
the lowest energy becomes the node 3. The average cur- 
rent during this time interval can be computed by assum- 
ing that a particle goes deterministically from the node 1 
to the node 3 via the path indicated by orange arrows in 
Fig. 3(b). After this, the protocol proceeds with a time 
interval during which a particle stays at the node 3 until 
the barrier configuration is such that the minimal span- 
ning tree is the one shown in Fig. 3(c). We then repeat 
our arguments until the protocol completes one period 
and the particle returns to the node 1. After the proto- 
col illustrated in Fig. 3(b)-(d), a unit current will be gen- 
erated through links {1,2}], {2,3}, {3,4}, {4,5}, {5,6}, 
{6,1}, and zero current will pass on average through links 
{1,4} and {2,4}. 
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VII. AVERAGE CURRENTS AND RATIONAL 
QUANTIZATION 

In this section, wc study the low-tcniperaturc adiabatic 
limit of Q(s) in the situation with permanent degeneracy 
(i.e., rational quantization setting). This is achieved by 
applying the low temperature — > oo limit to the adia- 
batic limit Q{s; (3) of the average current [Eqs. (31) and 
(32)], and ultimately to Eq. (54). We will further derive 
the averaging formula [Eq. (21)] that expresses Q{s) as a 
result of equal- weight averaging of a properly chosen set 
of integer currents. 



A. Minimal spanning trees in the degenerate case 

The minimal spanning tree, Xw, associated only with 
non-degenerate barrier data W , as established in sub- 
section VIA and appendix A maximizes its Boltzmann 
weight, and in the low-temperature /3 ^ oo limit domi- 
nates it. 

In the case of degeneracy of two or several barriers, 
we can instead identify a set of spanning trees that play 
a similar role to the role that minimal spanning trees 
play for graphs without barrier degeneracies. We will 
still refer to them as minimal spanning trees, having in 
mind that in the degenerate situation under discussion a 
minimal spanning tree is generically not unique. 

A minimal spanning tree can be obtained from the orig- 
inal graph by the following procediire. We remove links 
with the highest barriers one-by-one, and stop when the 
removal of any additional link of the chosen group breaks 
the graph into disjoint components. We then repeat the 
procedure with set of next-highest links and so on until 
wc end up with a tree. In a generic situation a minimal 
spanning tree is not unique, since its construction explic- 
itly depends on a particular choice of the withdrawn links 
within any group of degenerate links. 

For a non-degenerate (in the sense of rational quan- 
tization setting, i.e., no degeneracies in addition to per- 
manent degeneracies) barrier data W the set of minimal 
spanning trees is denoted by X-w By the construction of 
the minimal spanning trees, the set Xw depends on the 
ordering of the barrier groups heights, rather than their 
actual values. Therefore, similar to the integer quanti- 
zation setting, a notation Afj^j is in place. We will use 

a notation nwr{X) = csLTd{Xw) for the number of min- 
imal spanning trees of the graph X with given barrier 
data W. 



B. Properties of minimal spanning trees and 
associated graphical structures 

The crucial role of the sets Xw of minimal span- 
ning trees for the low-temperature limit is determined 
by their following properties: (i) The energies W{X) of 



all X £ Xw are identical, and (ii) the energy W{Y) 
of any spanning tree Y that does not belong to X\v 
is strictly larger then the energy of the spanning trees, 
and, therefore, the minimal spanning trees dominate the 
Boltzmann distribution in the low-temperature limit. 

All minimal spanning trees can be described using the 
following simple construction. For given barrier data 
W, let us order the barrier values M^aii > . • . > Wa„ 
and withdraw the links a G 2ti with the highest barri- 
ers Wa = Wst^. The resulting graph is generally dis- 
connected and is represented by a disjoint union of its 
connected components, referred to as children of X, with 
X naturally treated as their parent. We apply the above 
procedure to the children repeatedly, terminating when 
a child graph is a tree. All subgraphs obtained on the 
described way, including X are referred to as the derived 
subgraphs of X. Let Gx{W) be a tree graph whose nodes 
are the derived subgraphs of X . and the links represent 
the parent-child relation. Obviously the root of Gx{W) 
is given by X, whereas all leaves are represented by trees. 

With any non-tree derived subgraph we can associate 
its reduced counterpart, referred to as a reduced derived 
subgraph by contracting to a. point all links cixcept for the 
ones with the maximal barrier value. Naturally the nodes 
of a reduced subgraph are represented by the children of 
the original subgraph, with the links given by the links 
of the original subgraph with the maximal barrier value. 
The reduced subgraph of a tree [i.e. a leaf of ^x(W')] 
is defined as a result of contracting the tree to a point. 
Let Gx(W) be the tree obtained by replacing the nodes 
of GxiW), represented by derived subgraphs with their 
reduced counterparts. Note that all links of any reduced 
derived subgraph have the same barrier values, and the 
set of links Xi of our original network X is given by 
a disjoint union of the link sets of the reduced derived 
subgraphs and the links of the leaves of GxiW) (i.e. tree- 
like derived subgraphs). 

Naturally, any minimal spanning tree X G X^ is to- 
tally determined by a particular choice of a spanning tree 
for any reduced derived subgraph of X. The links of a 
minimal spanning tree are represented by a disjoint union 
of the link sets of the chosen spanning trees in the reduced 
derived graphs, and the link sets of the tree-like derived 
subgraphs [leaves of GxiW)]. 

Property (i) of minimal spanning trees follows from 
the general fact that all spanning trees of a graph have 
the same number of links, applied to the chosen spanning 
trees of the reduced derived subgraphs of X. Therefore, 
any minimal spanning tree has the same number of links 
of any given barrier value, which implies the same value 
of W{X) for all X G Xw- Property (ii) is validated in 
appendix A. 

The explicit description of the minimal spanning trees 

presented in this subsection also provides a convenient 
expression for the number nwi^) of the minimal span- 
ning trees as the product of the spanning trees numbers 
over all reduced derived subgraphs of X [i.e., nodes of 
Gx{W)]. 
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We conclude this subsection with noting that the pre- 
sented construction based on reduced derived subgraphs 
has a transparent physical meaning. Since in the (3 oo 
limit population relaxation within child subgraphs occurs 
much faster than the population exchange between them, 
the rcxiuced derived subgraphs represent the hierarchy of 
the relaxation processes time scales. 



the integer case, we obtain 

0(5) = E 



^ ^ ^ Q.fe'(^),(60) 



and, combined with the definition of Qij{X), Eq. (60) 
yields rational- valued Q{s). 



C. Low-temperature limit for the current and 
minimal spanning trees 



Now we show that in the case with permanent degener- 
acy, for any good protocol s, i.e. when s(t) stays within 
A4x at all times, there is a well-defined low-temperature 
adiabatic limit Q{s) = limi^^^limT-jy^oo Q{s', /3,td), so 
that Q{s) is rational-valued, and present explicit ways to 
evaluate this limit. Our derivation is very similar to the 
one for the integer case, presented in subsection VI B, and 
is based on identification of the low temperature limit of 
Eq. (54), where the differences refiect just the differences 
in the low-temperature limiting behavior of the node and 
tree Boltzmann distribution p^{x;f3) and q^{x:[3) for 
X & Uq (no energy data degeneracy) and x £ Ui (no bar- 
rier data degeneracy), respectively. The latter have the 
form 



hm p^{E;p)=u^' Yl ^' 
keiCE 



(57) 



for the node Boltzmann distribution, where ICe € 
denotes the subset of nodes k with the lowest value of 
El-, and Ue = card(/C£;) is the corresponding number of 
nodes, and 



lim Q''{W;p)=n^ V X 



(58) 



for the tree Boltzmann distribution. 

Following our derivation of the integer case we parti- 
tion the time period of the driving protocol into a set 
of alternating 0- and 1-segments exactly as described in 
subsection VI B and apply the low-temperature limit to 
Eq. (54). By exactly the same reasons as in the integer 
case, the integral terms vanish in the /? — )• oo limit. We 
further note that, since the ends s(to),s(t^) e Uq (1 U\ 
for all a, the low-temperature limits [Eqs. (57) and (58)] 
apply to the contact terms in Eq. (54). Therefore, denot- 



ing for brevity Ka = K.e(t^), IC'^ 
and u'^ = i'E(r'j we arrive at 



K 



""Y. X®{k'-k). (59) 

Finally, applying the bilinear form Q with the compo- 
nents Q''-^ to Eq. (59), exactly in the same way as for 



D. Averaging formula for rational quantized 
currents 

The expression [Eq. (60)] for the low-temperature adia- 
batic limit of the pumped current, not only demonstrates 
the rational nature of the pumped currents, but also pro- 
vides an algorithm for its computation. Yet, the afore- 
mentioned expression represents the pumped current as 
a rational superposition of integer currents generated at 
the 1-segments, so that the decomposition can depend on 
a specific choice of partitioning the time period into the 
set of alternating 0- and 1-segments. In this subsection 
we derive an expression for Q{s) in the form of a rational 
superposition of integer-valued currents, defined globally, 
i.e., for the complete period of the driving protocol. 

Each integer-valued component in the superposition is 
associated with a global degeneracy resolution. To intro- 
duce a degeneracy resolution we denote by [E] and [W] 
for any non-degenerate energy E and barrier W data 
the connected components of E and W \n Uq C Mx 
and U\ C M.X, respectively. Stated differently, [E] is 
fully determined by the lowest-energy subset /Ce, and 
we can consider [E] e X^; naturally \W] is fully deter- 
mined by a linear ordering of Xi, so that we can consider 
\W] € 0{Xi), with 0{Xi) being the set of permutations 
(or equivalently linear orderings) of ^"1. A global degen- 
eracy resolution (k, X) associates a node kE G )Ce with 
any [E] e Xq and an optimal spanning tree Xw € Xw 
with any [W] € 0{Xi). Obviously, the number nx of 
possible degeneracy resolutions is given by 



nx = Yl I'E JJ nw 

[E](^Xo [W]&0(Xi) 



(61) 



To derive the AF, we consider Eq (60) and we add the 
identity resolution, 



v^nx 



EE E E 1 = 1. 



(62) 



to each term in the sum over time intervals a. Using 
that 'Y^]^E]^Xo 1 = ^-iid that the sum over all possi- 
ble identity resolutions, given by '}2(k x) = 



j[E]eXo l^[w]eO(Xi) 



we obtain 



Q{s) 



= E^ Eo.„.i(^). 



(63) 



(fe,X) 



Since ^ is a constant, we can interchange positions of 
the sums in (63). We then collect the piecewise currents 
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Qkak' (X) in groups. The sum of the currents in each 
group represents some integer valued conserved (in fact, 
cyclic) current that starts and ends in the same node. 
Each such a distinct integer-valued conserved current is 
associated with a corresponding global degeneracy reso- 
lution (fe,X), so that 

Qis; k,X) = J2 i^wj, (64) 

a 

where the summation runs over the 1-segments for some 
proper partition of the time period into a set of alter- 
nating segments, and the initial and final nodes coincide. 
Although the definition of the integer reference current 
Q{s; k, X) involves a partition into alternating segments, 
we will argue in the following subsection VIll E that it 
does not actually depend on a particular partition choice, 
and therefore is fully determined by a degeneracy resolu- 
tion and the driving protocol. 

From Eqs. (63) and (64) we finally arrive at the aver- 
aging formula for the pumped current 

^ (fe,x) 

that expresses Q{s) as a result of averaging the integer 
cyclic currents over all degeneracy resolutions with iden- 
tical weights. 

The fact that the integer currents in (65) are cyclic 
implies that the total pump current, Q{s), is conserved, 
as was expected. 

VIII. WINDING-INDEX FORMULAE FOR 
QUANTIZED CURRENTS 

In this section we present a winding index representa- 
tion of the pumped current that provides an alternative 
way of calculating the low-temperature adiabatic limit 

A. Winding index and intersection index 

The quantized currents in both integer and rational 
cases can be represented in terms of the winding numbers. 
As it was briefly discussed in subsection 11 D, the space 
Aix — M.X \ Dx can be viewed as a result of withdraw- 
ing a set Dx of bad parameters, the latter referred to as 
the discriminant set. The discriminant set Dx consists 
of points {E, W) e Dx with simultaneous degeneracy of 
some energies Ej = Eji and barriers W% = Wa', and 
therefore has codimension codmi[Dx) = 2. Therefore, 
for a good driving protocol s that by definition avoids 
Dx, one can think in terms of winding, i.e. how many 
times the driving protocol winds around the discriminant 
set, such a number is usually referred to as the winding 
number. To formalize the notion of the winding num- 
ber, one spans a surface A on s, i.e, we have for its 



boundary dA = s, and coimt the intersections of A with 
Dx, weighted with the proper sign factors, depending on 
the relative orientation of A and Dx at the intersection 
point, the obtained integer number, denoted A ■ Dx be- 
ing known as the intersection index. The winding index 
is defined by s * Dx = A ■ Dx ■ 

There are two things that should be taken care of to 
build a winding-index representation of the generated 
current. First of all, the winding index shoiild not de- 
pend on a specific choice of the surface A spanned on 
s. This will be achieved if Dx is a cycle, i.e., it does 
not have a non-trivial boundary. Second, the winding in- 
dex should be represented by a conserved current on our 
graph X with rational components, rather than by just a 
single integer number. This will be achieved by extend- 
ing Dx to a cycle Vx with coefficients in the space of 
conserved currents with rational components, so that we 
have Q{s) = s^Vx = A-Vx- The cycle Vx is naturally 
referred to as the current-generating cycle. 



B. Geometrical structure of the discriminant and 
the current-generating cycle 

The discriminant set Dx has a very natural decom- 
position into cells, so that with each cell we can asso- 
ciate an elementary conserved current on our graph X. 
The set of elementary currents, associated with the cells 
totally determine the current-generating cycle 2?x- An 
elementary cell is labeled by the following data: an un- 
ordered pair {J, J'}, related to a lowest-energy degen- 
eracy Ej = Eji, and an unordered pair {21,21'}, related 
to a barrier degeneracy W% = W%i , together with an el- 
ement p G 0(A'i/{2l,2l'}), i.e., an ordering <Bi < . . . < 
Q3j < {21, 2t'} < 2lj+2 < . . . < *B„ "with degeneracy" 
of Xi = {*Bi, . . . , «Bj, 2t, 21', 2lj+2, • • • , *Bm}. An open cell 
Djj'MW,p consists of the parameter sets {E,W) that 
satisfy the following constraints: Ej = E'j < Ex for any 
I ^ J, J', and W<8, < ... < W2}, < W^a = W<^, < 
W(Sj^2 < ■ • • < W<s,„. Stated differently in an open cell 
exactly two specified lowest energies are degenerate, and 
exactly two specified barriers are degenerate, and the or- 
dering of all barriers is also fixed. Any open cell is con- 
tractible and topologically equivalent to an open multidi- 
mensional disc (ball). A closed cell ast'.p, obtained 
as a topological closure of the corresponding open cell, is 
defined by the modified conditions, where the signs "less" 
are replaced with "less or equal". Naturally, different 
open cells do not intersect, whereas Dx = U-^J'J'',2ia',p- 
Stated differently, the discriminant Dx can be viewed as 
a bunch of closed cells, glued together along their bound- 
aries. 

There is a very natural and intuitive way to associate 

with each cell D j j-'^aa'.p a- conserved current V j j\%'Si\p 
by choosing a reference point {E,W) € Djj/^<^<2i>^p, 
considering a good driving protocol s that winds once 
around the cell staying close to the reference point, fol- 
lowed by evaluating the generated current and setting 
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T^jj'.%^',p = Q('S)- More specifically, we choose s{t) 
by setting all parameters, except for Ej and Wa being 
time-independent, whereas 

Ej = Ej, + £cos(27rT), Wa = Wa' + esin(27rr)(66) 

with £ being small enough not to cause any additional 

degeneracies. The generated current Q{s) can be readily 
calculated, by partitioning the time period [—1/8,7/8] 
into four alternating 0- and 1-segments [—1/8,1/8], 
[1/8,3/8], [3/8,5/8], and [5/8,7/8]. This results in 



V 



1 



E E 



X g {Qkk'{X) + Qk'k{X')), (67) 

X,X' 

where the X, X' summation goes over X e <^w„,„,, and 

X' e Xw„,,„, . The set of conserved rational- valued cur- 
rents T>jji_%<^i,p on X, associated with the open cells 
Djj'^f^fU'^p, will be referred to as the current-generating 
cycle, and denoted Vx- It is convenient to represent the 
current-generating cycle as a formal sum 



V 



X 



(68) 



over the open cells lu = {J^J^',^'Ql',p) that decompose 
the discriminant set Dx- 

The derived expression [Eq. (67)] for the elementary 
currents has a simple intuitive interpretation. For a 

and 



aa'.p' 



given choice of k € J , k' € J' , X G Xw, 

X' € '^w'si'si p concatenate the paths lkk'{X) and 
lk'k{X'). This results in a closed path in X that gener- 
ates a conserved integer- valued current. Averaging these 
currents over all allowed choices of {k,k' ,X ,X') with 
identical weights results in the elementary current V,^, 
associated with the cell w. 



C. Pumped current as winding index 

To define the winding index s * Vx with values in con- 
served currents on X having rational components, wc 
span the surface ^ on s (more specifically a disc), so that 
it avoids the boundaries D^^ — of the cells D^^, i.e., 
the intersections arc restricted to the open cells and 
are transversal, in particular, the intersection A fl Dx is 
represented by a finite number of isolated points. Using 
the representation of Eq. (68) and the general concept of 
winding, briefly described in subsection VIII A, we define 

s*Vx = A-Vx = Y,'^^(-^-^^)^ (69) 

where u = {JJ',^^',p), and the standard, i.e., integer- 
valued, intersection index is defined by 



A- D^= ^ inda 



(70) 



and the intersection index indj, = ±1, associated with 
an intersection point, reflects the relative orientation of 
A and at the intersection x. Since a change in ori- 
entation of is accompanied by a sign change of the 
elementary current V^^, the definition given by Eq. (69) is 
invariant with respect to re-orientations of the cells D^^. 

The winding-index formula for the low-temperature 
adiabatic limit of the pumped current 



Q{s) = s * T>x 



(71) 



provides an alternative way to calculate Q{s). One can, 
e.g., choose some reference point y G A^x and consider 
A as the cone of the driving protocol s with the origin in 
y. The obtained surface A spanned on s is smooth every- 
where except for the cone origin y, which is not a problem 
since by construction y is off the discriminant set Dx and 
does not participate in the intersection. In a generic situ- 
ation the intersection A H Dx is transversal, however, to 
convert the presented idea into an algorithm, one should 
introduce a specific procedure of resolving non-transverse 
intersections. This issue will not be addressed here. 

The winding-index expression [Eq. (71)] for Q{s) is 
very intuitive and can be rationalized as follows. Due 
to the robustness property of Q{s) and given a spanned 
surface A with a transversal intersection A n Dx^ the 
driving protocol s can be replaced by a set Sx of driv- 
ing protocols that wind around the intersection points 
X € Af\ Dx, so that the pumped current remains the 
same, i.e., Q{s) = Y^xeAciDx Qi^^)- ^J' observing that 
Qi^x) = T^ui^, where uj^ labels the cell that contains 
X e D^^ and making use of Eq. (69), we arrive at the 
winding index formula [Eq. (71)]. 



D. Cyclic nature of the current generating cycle 
and the space of robust parameters 

As stated above for the winding index s * Vx to be 
defined properly, i.e., for the intersection index A ■ Vx 
to be independent of a particular choice of the surface A 
spanned on s, the object Vx should be a cycle, i.e., it 
should not have a non-trivial boundary in some proper 
sense that has not been identified yet. On the other 
hand, the winding-index representation [Eq. (71)] clearly 
demonstrates the invariant nature of s * Vx, i.e., the in- 
dependence of the choice of surface A, which means that 
referring to Vx as a cycle with coefficients in conserved 
currents on X definitely makes sense. 

To put the above argument on a more formal basis 
we note that, since our cells D^^ have codimension 2, 
their boundaries have codimension 3, and there are three 
types of boundaries: (i) triple energy degeneracy together 
with a simple barrier degeneracy, i.e., Ej — Eji — Ejn 
and W% = Wai', (ii) simple energy degeneracy together 
with a triple barrier degeneracy, i.e., Ej = Ej' and 
Wgi = Wa' = Wa", and (iii) simple energy degener- 
acy together with a pair of distinct barrier degenera- 
cies, i.e., Ej = Eji, Wa = M^a', and = W^i. Of 
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course, wc need to specify the orderings of the degenerate 
groups of barriers together with the non-degenerate bar- 
riers, which represents the total boundary as a union of 
closed cells of codimension 3, whose intersections corre- 
spond to higher-order degeneracies, having codimension 
4, and will not play role in our consideration. A type- 
(i) boundary cell has 6 current-carrying cells attached to 
it, e.g., Ej = Ej' < Ejii and W% = W%', with the 
ordering with respect to the other barriers being deter- 
mined by their ordering in the considered boundary cell. 
The situation can be illustrated in a 3-dimensional space 
of relevant parameters (Eji — Ej). [Ejn — Ej), and 
(M^ai' — W%). Similarly, a type-(ii) boimdary cell has 6 
current-carrying cells attached to it. The relevant pa- 
rameters are {Ej> — Ej), {W<^i — VFa), and (PFa" — W'a). 
A type-(iii) boundary cell shares the boundary of 4 cells, 
e.g., Ej = Eji, T'Vgi < VFsi', and W<s = VT®', with 
the relevant parameters {Eji — Ej), {W%' — Wa), and 
(W2}' -W<s). 

The independence of the weighted winding index s * 
Vx = A ■ Vx on a particular choice of the spanned sur- 
face, which actually takes place due to Eq. (71), boils 
down to the condition that for each boundary cell, the 
sum of elementary currents I?^^, associated with the cells 
D^^ attached to it, should be equal to zero. Of course, all 
should be accounted for with the proper signs that, 
e.g., correspond to all elementary currents to be incoming 
to the boundary. Stated differently we can introduce the 
boundary operator that associates with Vx = Yliu ^ 
formal sum over the boundary cells Bg, labeled by q: 

dVx = dJ^'D.oj = J2^^9^ = J2eJ2'^-^ (^2) 

WW a ujDq 

where uj D g naturally means 5^ D Bg, and Eq. (72) is 

consistent with a natural definition du! = X^pcw ^- Using 
this notation we can reformulate the requirement of the 
winding index s * Vx to be well-defined in two equivalent 
forms (invariant and component) 

&Dx=0, J2V^=0, (73) 

with the invariant form claiming that the boundary of 
Vx is trivial, i.e., Vx is a cycle, and the component 
form to be some kind of a higher dimensional Kirchhoff 
law, claiming that the sum of currents, represented by 
V^ that are incoming to any boundary component Bg is 
equal to zero. 

The derived winding index representation allows the 
subspace Aix of good parameters to be readily extended 
to a space Mx ^ -^x, described as the subset of ro- 
bust parameters Mx C Mx, i-e, that provides robust 
behavior of the low-temperature adiabatic current. Not- 
ing that only those cells D^^ with V^^ ^ participate in 
the current generation, we have: 

M'x=Mx-D^x, D^x= \J D^- (74) 



E. Topological view and invariant nature of the 
reference currents 

In this subsection we present some simple and intuitive 
topological arguments that support the assertion that a 
reference current Q(s; fc, X) is independent of the par- 
ticular choice of the partition of the protocol period into 
alternating segments. We will also apply them to demon- 
strate that the pumped current Q{s) is robust, i.e., it does 
not change upon continuous deformations of the driving 
protocol s, as long as s entirely stays in the subspace M.x 
of good parameters. Our arguments are based on intro- 
ducing an alternative invariant definition of the reference 
currents Q(s; fc, X) that does not involve any partitions. 

To that end consider the space M.x x l-'^l, where \X\ 
denotes the geometrical realization of our graph X, i.e., 
the nodes of X are represented by points in |A"|, whereas 
the links are represented by segments, whose ends are at- 
tached to the points that represent the boundary nodes. 
In simple words |X| is a one-dimensional space, obtained 
when we "draw" a graph. We can consider a stochas- 
tic trajectory of a particle as a trajectory of a point 
(a;, Tj) e Mx^ \X\, where the component x moves in the 
"good" parameter space Adx deterministically according 
to the driving protocol, whereas rj moves in a stochastic 
manner in the graph space \X\. We fix a degeneracy 
resolution (fc, X) and further impose two restrictions on 
the particle motion: (i) when x = {E, W) € Ui, i.e., the 
barrier data is non-degenerate, the particle position rj on 
the graph is restricted to ij E \Xw\ C l^l, (ii) when x 
docs not belong to Ui, i.e., the barrier data degeneracy 
occurs, and hence x G Uq, the particle should be located 
at the node fc^;, i.e. r] = ks- 

Consider a periodic stochastic trajectory (s(T),J7(r)) 
in the restricted space, where s(r) is the driving pro- 
tocol over a time period. The r/(r) component of the 
trajectory can be arbitrary, provide the above restric- 
tions are satisfied. Despite a large variety of allowed 
trajectories, the currents generated by all these trajec- 
tories are identical. There is a topological reason for 
that: the particle motion along the graph X is always re- 
stricted to a tree: J7(t) € when s(t) G Ui, and 
rj{T) = kEir), otherwise. Since trees are contractible, all 
trajectories J7(t) for the same driving protocol are topo- 
logically (homotopically) equivalent, and hence produce 
the same current, denoted Q{s;k,X). To demonstrate 
that the current Q{s; k, X), defined above in an invari- 
ant manner using simple topological terms, is identical to 
the one defined by Eq. (64) in terms of a partitioning, we 
evaluate the "topological" current for a special choice of 
the trajectory r/, so that {s,ri) goes along the restricted 
subspace. For a 0-segment we set ■r]{T) — ksir), for a 1- 
segment la the trajectory r7(r) connects ksiTa) to kEiT^) 
via the path lkE{Ta)kE{T'){Xwa)- By inspecting the def- 
initions we see that the "topological" current generated 
by the described trajectory r; reproduces Eq. (64). 

The topological view of the reference currents 
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Q(s; k, X) combined with the averaging formula allows 
the robustness property to be demonstrated in a simple 
and intuitive way. Indecid, according to the above a ref- 
erence current Q{s; k, X) can be obtained by lifting the 
closed path s to the restricted subspace of Mx x \X\ 
(the lift being unique up to a homotopy), followed by 
projecting the lift into |X| with evaluating the current, 
associated with the obtained path in the graph X. All 
operations preserve homotopy: if s and s' are homotopic, 
i.e., s can be deformed to s' in a continuous manner, then 
so are their lifts, and further the projections of the lifts 
into \X\. Since the conserved current associated with a 
closed path in the graph X depends on the homotopy 
equivalence class of the path only, any reference current 
is a homotopy invariant of the driving protocol s, i.e., 
is robust. This implies the robustness of Q{s) via the 
averaging formula [Eq. (65)]. 



F. Winding index representation and current tubes 

Spanning a surface A on the driving protocol s and 
applying the Stokes theorem, we can express the current 
Q{s;(3) generated at finite temperature 

Q{s;l3) = / dAhdp^ = [ F 
J A J A 

= V / F"^ {E, W; (3)dEj A dWc,, (75) 

as a surface integral of the (antisymmetric) curvature ten- 
sor F that takes values in the vector space of conserved 
currents on X and has the form 



The representation of Eq. (75) expresses the adiabatic 
(and not necessarily low-temperature) current as a flux of 
F through a spanned surface A, and therefore provides a 
physical interpretation of the winding index representa- 
tion of the low-temperature limit. Indeed, Eq. (71) with 
the winding index defined by Eq. (69) as the weighted 
intersection index A ■ Vx , the low-temperature current 
can be viewed as the flux of some curvature F, which 
is restricted to the discriminant set i?^, and hence is 
described by a (5-function type distribution. This is not 
surprising, since it describes the low-temperature /3 — >■ oo 
behavior of the curvature tensor F, whose explicit ex- 
pression for arbitrary /3 is given by Eq. (76). For low 
temperature, F is substantially concentrated in a nar- 
row tube around , with the tube transverse size tend- 
ing to zero exponentially with /? — > oo, whereas the flux 
through a small transverse surface that intersects at 
a point X e reaching a finite limiting value equal 
to r>u)- This can be readily established by comparing 
Eq. (76) with the the explicit expressions [Eqs. (31) and 
(32)], provided by the Kirchhoff theorem, which results 



F°'\E,W;P) 



k,X 



kX 



idQ%{W,P)dpB^E,p) 



dEi 



(77) 



followed by a straightforward study of the low- 
temperature behavior of the weighted curvature F. 

Naturally the cyclic property dDx = of the cur- 
rent generating cycle Vx is nothing else than the 
low-temperature counterpart of corresponding property 
dF{(3) — 0, the latter holding for all temperatures. 



IX. DISCUSSION 

We have considered the current generated in a stochas- 
tic network due to periodic time-dependence of the ki- 
netic rates. We demonstrated that for a typical driving 
protocol its low-temperature adiabatic limit Q{s) is well- 
defined and robust, i.e., does not change at all with re- 
spect to moderate perturbations of the driving protocol. 
We established that the integer quantization of this cur- 
rent follows from the fact that in the considered low tem- 
perature limit no average current flows over the largest 
barrier if the corresponding link can be removed without 
breaking graph into disjoined component. This property 
holds true even if an alternative path forces the system to 
overcome a series of barriers and even though adiabatic 
limit assumes that system performs many random tran- 
sitions through all graph links. As a result, effectively, 
motion on the graph is restricted to a set of spanning 
trees with identical weights. Using this observation, we 
derived an explicit expression for the value of the frac- 
tional quantized current in a general graph and for a 
general protocol. 

Our main result - the averaging formula, given by 
Eq. (65), provides a combinatorial approach to calcula- 
tion of the fractionally quantized currents. Interestingly, 
it tells that fractional quantization is reduced to integer 
quantization case, namely, the final result corresponds 
to averaging over all properly defined legitimate degen- 
eracy resolutions of the barriers and node energies, with 
the individual contributions, labeled by the degeneracy 
resolutions being represented by integer quantized cyclic 
currents that can be found using a deterministic picture 
for integer current calculation. This simplification can be 
used, e.g., for efficient numerical calculations of fraction- 
ally quantized currents in large networks, avoiding time- 
expensive matrix inversion and diagonalization that are 
used in pump current calculations [29, 35, 36]. 

We have also established the winding index view of the 
current quantization in the low-temperature adiabatic 
limit. This is achieved by identifying the cellular struc- 
ture of the space D^-^ C Mx of non-robust parameters, 
as well as associating with each cell of maximal dimen- 
sion (which corresponds to codimension 2) a conserved 
current, so that for each cell of codimension 3 the sum of 
the currents (with a proper sign choice) , associated with 
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the codimension 2 cells attached to it, is zero. The above 
identification leads to a notion of the current generating 
cycle Vx with the coefficients in the conserved current, 
so that the current Q{s) generated by a robust driving 
protocol is given by the winding index of s and Vx, as 
given by Eq. (71). The winding index representation is 
reminiscent of the Aharonov-Bohm effect in quantum me- 
chanics, or strictly speaking its generalization to the case 
of arbitrary dimension of the space Aix, rather than just 
3. In the latter case the space Dx is represented by a 
graph, with the conserving currents residing on the links 
and the current conservation at the nodes (the sum of 
currents entering any node is equal to zero). We showed 
how to determine the space Dx explicitly, and provided 
a procedure that assigns a quantized conserving current 
with each cell of D^, resulting in the cycle Vx- Our 
work provides a direct algorithmic approach for a com- 
plete classification of fractionally quantized currents in 
any specific model of stochastic Markovian motion of a 
mesoscopic system under an influence of a periodic driv- 
ing protocol. 

We conclude with noting that in section IV we imple- 
mented the KirchhofF theorem to develop a useful tool 
for applying the low-temperature adiabatic limit, i.e., we 
derived an explicit expression [Eqs. (27) and (28)] for the 
average generated current for a general situation, i.e., 
finite values of the inverse temperature /3 and driving 
protocol period td. This expression might have its own 
value, especially in the adiabatic td ^ oo, still finite- 
temperature case, where the expressions [Eqs. (31) and 
(32)] show some surprising symmetry between the sets 
of nodes and spanning trees of our original network X. 
This symmetry could be a manifestation of some hidden 
duality, with the nodes of the dual graph defined as the 
spanning trees of the original one. We note that on cyclic 
graphs such a duality is well known [37]. In particular, 
it was used to derive unusual fluctuation relations for 
currents [25]. The existence of such duality in models 
defined on generic graphs remains an open question and 
could be subject of future work. 
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Appendix A: Properties of minimal spanning trees 

In this appendix, we rationalize the properties of the 
minimal spanning trees, described in subsection VI A for 
the integer quantization setting, as well as in subsec- 
tion VII B for rational quantization setting. More specifi- 
cally we are dealing with property (i) in the integer case, 
and property (ii) for the rational quantization setting. 



We reiterate that the difference between the two cases 
is that in the integer case the minimal spanning tree is 
unique, and is denoted by X^, whereas in the rational 
case there is a set of minimal spanning trees, denoted 
by Xw with all minimal spanning trees having the same 
value of the energy function W{X). The property un- 
der consideration can be formulated in a unique fashion 
for both integer (non-degenerate) and rational (degener- 
ate) cases: The value of the energy function of a minimal 
spanning tree is strictly lower than the energy value of 
any spanning tree that is not minimal. Although the 
non-degenerate (integer) situation is a particular exam- 
ple of the degenerate (rational) counterpart, for the sake 
of presentation clarity we start with the non-degenerate 
situation, followed by demonstrating how the presented 
arguments can be extended to the degenerate case. Note 
that in both cases it is sufficient to show that if a span- 
ning tree X is not a minimal spanning tree, then there is 
a spanning tree X' with W{X') < W{X). 

We start with the non-degenerate case. Let 
{ai, . . . ,ak) and (/3i, . . . , /3fc) be the elements of Xi \ 
{Xw)i and Xi \ Xi, i.e., the links withdrawn to ob- 
tain the trees Xw and X, respectively, in the decreasing 
order of their barrier levels. The latter means Wai > 
. . . > Wock and > . . . > Wp^^ . Let j be the smallest 
index such that aj ^ Pj, and ai = Pi for i < j. Con- 
sider a spanning subgraph X' C X, obtained from X 
by withdrawing the edges ai,...,aj_i, or equivalently 
. . . , Pj-i- By the definition of X\y the edge aj is not 
a bridge of X' (i.e., its withdrawal does not disconnect 
the graph), and W/s^ < Wa- for i > j. Let and X'^ 
be the two trees obtained from X by withdrawing the 
edge aj. Then there is at least one edge, say [Ss among 
/3j, ■ . . ,(3k that connects X^ to X^, since otherwise the 
edge aj would be a bridge of X'. Therefore, replacing the 
edge f3s with aj in X results in another spanning tree, 
denoted X' with W{X') < W{X), since VF„. > Wp^. 

The above arguments may be extended to the degener- 
ate case in a straightforward way. To that end, we note 
that since the number of links in each group S2t C Xi 
of degenerate barriers that are withdrawn to obtain a 
minimal spanning tree is the same for all minimal span- 
ning trees, we can consider the groups that provided 
at least one link to the withdrawal procedure and put 
them in the decreasing order (2li, . . . , 2lfe) of barrier val- 
ues > . . . > Waifc • We also denote by (m, . . . , n^) 
the numbers of links withdrawn in each group to obtain 
a minimal spanning tree (as noted above, these num- 
bers are the same for all minimal spanning trees). We 
further consider a spanning tree X that is not minimal 
(i.e., does not belong to Xw) and denote by (*Bi, . . . , 23;) 
and {rui, . . . ,mi) the groups (in the decreasing order 
of the barrier levels) and the numbers of links, respec- 
tively, that provided to the withdrawing procedure that 
resulted in X. Let j be the smallest number so that 
(2lj,nj) 7^ (*Bj,mj), and (21^,71,^) = {^i,m,i) for i < j. 
Let X' C X he a. spanning subgraph, obtained from X 
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by withdrawing all edges that belong to |Ji=i \ -^i- 
Since, by the construction of minimal spanning trees, 
card(2lj \Xi) < rij, there is always a link a e 2tj, which 
does not bdong to X[ \ Xi, and therefore Wa > Wp 
for all (3 e X[\Xi. Similar to the integer case we can 



replace the link a that belongs to Xi with some prop- 
erly chosen link j3 from X[ \ Xi to obtain a tree X' with 
W{X') < W{X). 
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